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Abstract 

The design and generation of modified involute helical gears that have a localized and 
stable bearing contact, and reduced noise and vibration characteristics are described. The 
localization of the bearing contact is achieved by the mismatch of the two generating surfaces 
that are used for generation of the pinion and the gear. The reduction of noise and vibration 
will be achieved by application of a parabolic function of transmission errors that is able to 
absorb the almost linear function of transmission errors caused by gear misalignment. The 
meshing and contact of misaligned gear drives can be analyzed by application of computer 
programs that have been developed. The computations confirmed the effectiveness of the 
proposed modification of the gear geometry. A numerical example that illustrates the devel- 
oped theory is provided. 



Nomenclature 


a Parabola parameter (fig. 2(a)) 

a c Modification coefficient of pinion rack-cutter surface (fig. 5) 
b Slope of linear function (fig. 1(b)) 

E pc The shortest distance between the pinion and rack-cutter E e (fig. 6(b)) 
E pg The shortest distance between the pinion-gear axes (fig. 7) 

E gt The shortest distance between the gear and rack-cutter E t (fig. 6(a)) 

m 2 1 Gear ratio 

M,j Coordinate transformation matrix (from Sj to Si) 
n r Unit normal vector to rack-cutter surface £ r (r = c, t) 

Unit normal vector to surface E, represented in coordinate 
system 5/ (t = p,g) 

Ni Number of teeth of the pinion (i = l,p) or the gear (i =2 ,g) 

N r Normal vector to rack-cutter surface £ r (r = c, t ) 

p n Circular pitch in normal section (fig. 3) 

P n Diametral pitch in normal section 

7 \ Radius of the pitch circle of the pinion (or gear) (i = p, g) 

r, -,r* Position vector of surface E* 

r^ Position vector of surface £, represented in coordinate system 5/ 

s, s r Displacement of rack-cutter E r (r = c,t) (fig. 6) 

Si Coordinate system i 

m,6i Surface parameters of E, (i = p,g) 
u r ,l T Surface parameters of E r (r = c, t) 

Velocity of rack-cutter surface point (r = c, t) 
yi'i) Relative velocity of surface E, point with respect to surface £_, point 
a 0 Normal pressure angle (fig. 3) 

/3 o Helix angle on the pinion (gear) pitch cylinder (figs. 3 and 4) 

6 Elastic approach of pinion and gear tooth surfaces 
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A£ 

AA 0 

Aq 

A 7 x 

A 7v 

A<£ 2 ,Ai /> 2 

A <^2 

A 0 

S, 


Change of center distance 

Change of pinion lead angle on the pitch cylinder 
Displacement of contact point caused by misalignment 
Misalignment angle formed by crossed gear axes (fig. 8(a)) 
Misalignment angle formed by intersected gear axes (fig. 8(b)) 
Transmission error (fig. 2) 

Vector of the angle of compensating turn of gear 2 

Lead angle on pinion pitch cylinder 

Pinion (i = p) and gear (i = g ) tooth surfaces 


E r Rack-cutter surfaces (r = c,t) 

Rotation angle of gear i (i = 1,2 ,p,g) (figs. 2 and 7) 

1 Pg t Rotation angle of gear being in mesh with the rack-cutter St (fig. 6(a)) 

ib c Rotation angle of pinion being in mesh with the rack-cutter E c (fig. 6(b)) 
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1. Introduction 


Conventional helical involute gears are designed for transformation of rotation between 
parallel axes. Theoretically, the gear tooth surfaces are in line tangency at every instant, 
along a straight line that is a tangent to the helix on the gear base cylinder. However, the 
line contact of gear tooth surfaces can be realized only for an ideal gear drive. In reality, the 
crossing of axes of rotation (instead of being parallel) and errors of lead angle result in the 
so-called edge contact, as a specific instantaneous point contact caused by curve-to-surface 
tangency. Here, the contacting curve is the edge of the tooth surface of one of the mating 
gears and the contacting surface is the tooth surface of the other one. 

Trying to avoid the edge contact, the manufacturers of helical gears use various methods 
of crowning (deviation) of the theoretical gear tooth surfaces. However, the applied methods 
of crowning have not been complemented with the analysis of transmission errors caused 
by misalignment. Our investigation shows that improper crowning may avoid edge contact 
but cannot avoid the appearance of transmission errors of the shape shown in fig. 1. The 
function of such transmission errors is piecewise, almost linear, and has the frequency equal 
to the cycle of meshing of one pair of teeth. The above mentioned transmission errors cause 
high vibration and noise and therefore such transmission errors must be avoided. This can 
be achieved by application of computer numerically controlled (CNC) machines that have 
opened new perspectives for generation of gear tooth surfaces with improved topology. 

The intent of this paper is to describe a modified topology of low-noise involute helical 
gears that satisfies the following requirements: 

(1) The noise and vibration of helical gears are reduced substantially by application of a 
predesigned function of transmission errors of a parabolic type (fig. 2). Such a function can 
absorb (see below) an almost linear function of transmission errors shown in fig. 1. 

(2) The bearing contact is localized. Theoretically, the tooth surfaces are in tangency at 
every instant at a point instead of a line. The contact of gear tooth surfaces at every instant 
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is spread over an elliptical area due to elastic deformation of gear teeth. The dimension of 
the instantaneous contact ellipse can be controlled by choosing proper design parameters. 

( 3 ) The proposed gear tooth surfaces can be generated by two rack-cutters designed for 
generation of the pinion and gear, respectively. A nonlinear transmission function in the 
process for gear generation must be provided and this can be accomplished by application of 
the CNC machine. A linear transmission function is provided in the process for the pinion 
generation. 


2. Interaction of Parabolic and Linear Function of Transmission Errors 

The ideal gears transform rotation with constant gear ratio 77*21 = -rjr, and the ideal 

■‘*2 

transmission function is 


«(#>) = (i) 

where N x and iV 2 are the tooth numbers of the pinion and gear, respectively. 

However, the crossing of gear axes (instead of being parallel), intersection of these axes 
and errors of lead angle cause a transmission function <£2(^1) that is shown in fig. 1 (a). 
Our investigation (see sections 4 - 6 ) shows that the function of transmission errors caused by 
above mentioned errors of misalignment is a piecewise almost linear function of transmission 
errors A <^2(^1) with the frequency of a cycle of meshing for one pair of teeth (fig. 1 (b)). 
Here: 


Ni 

A <^2(^1) = ^(Pi) 

iV 2 


( 2 ) 


Transmission errors of this type cause a discontinuity of the transmission function and a 
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big jump of the angular velocity of the driven gear at transfer points (when one pair of teeth 
is changed to another one). Therefore, vibration and noise become inevitable. 

It was proven [1,2,4] that a predesigned parabolic function of transmission errors in- 
teracting with a linear function will become a parabolic function with the same parabola 
parameter. A parabolic function of transmission errors is much more preferable than a lin- 
ear function since the transmission function will be a continuous one, the jump of angular 
velocity of the driven gear and the stroke at the transfer point will be substantially reduced. 

Fig. 2(a) shows the sum of two functions of transmission errors 


A<A 2 (<M = AfPifa) + = b<f> 1 — o.<f>\ 


( 3 ) 


The first one, Ad> 2 l ^(<£i), is caused by misalignment. The second one, Ad> 2 2 *(<?!>i), is a 
predesigned parabolic function which exists even if misalignment does not appear. It is easy 
to verify that equation (3) represents in the new coordinate system (AV> 2 , V’l) the parabolic 
function (fig. 2(b)) that is designated as 


At/> 2 = — aij)\ 


( 4 ) 


The parabola parameter a in equations (3) and (4) is the same. Axes of coordinate 
system (AV> 2 > V’i) and (A<£ 2 , <j>\) are parallel but the origins are different. The coordinate 
transformation from (A</> 2 , <f>i) to (At/> 2 , xp x ) is represented with the following equations 


At /> 2 = A <^2 ~ 


4a’ 




( 5 ) 
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The difference between functions A<£ 2 (<£i) and Ai/» 2 (V , i) is the location of the couple of 
points ( A,B ) and the respective points (A*, B m ) (fig. 2(a)). The symmetrical location of 
( A,B ) is turned into the asymmetrical location of (A*,B‘). However, the interaction of 
several functions AV’ 2 (V’i) determined for several tooth surfaces being in mesh may provide 
a symmetrical parabolic function of transmission errors as shown in fig. 2(b). This can be 
achieved if the parabolic function will be predesigned in the area (fig. 2(a)) 

<h(B)-MA)> |^ + 2c (6) 


where c = — . Requirement (6), if observed, enables to provide a continuous function 

2 a _ ' ' 

A0 2 (t£x) for the range of — where Ni is the pinion tooth number. It will be shown below 
(see sections 5 and 6) that functions of transmission errors caused by angular errors (such 
as the crossing and intersection of the axes of rotation, error of the lead angles) are indeed 
piecewise linear functions, and the coefficient b can be determined knowing the angular error 
caused by misalignment and the design parameters of the gear drive. 


3. Surfaces of Rack-Cutters 

The imaginary process of generation of conjugate tooth surfaces is based on application 
of two rack-cutters that are provided respectively by a plane E t and a cylindrical surface E c 
that differs slightly from plane E t (see fig. 3). The rack-cutter surfaces E t and E c are rigidly 
connected each to other in the process of the imaginary generation, and they are in tangency 
along a straight line, ObZt (fig- 5). This line and the parallel axes of the gears form angle 
/3 0 , that is equal to the helix angle on the pinion (gear) pitch cylinder. The normal sections 
of the rack-cutters are shown in figs. 3 and 5. Rack-cutter surface E c generates the pinion 
tooth surface E p , and the rack-cutter surface E t generates the gear tooth surface E s . 
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Gear Rack-Cutters E t 

Using figs. 3, 4 and 5, we represent the transformation matrix from system S a to S T 
(r = c, t) and S b to S a as follows 


M ro = 


10 0 0 
0 cos 0 O sin 0o 0 
0 — sin 0 O cos 0 O 0 

0 0 0 1 


(7) 


A^ai = 


cos a 0 — sin a 0 0 —d p cos a 0 
sin cc 0 cos a 0 0 a m — d p sin a„ 
0 0 1 0 

0 0 0 1 


( 8 ) 


M rb = 


cos a 0 

sin a 0 cos 0 O 
— sin a 0 sin 0 O 


— sin a 0 0 

cos a 0 cos 0 O sin 0 O 

— cos a 0 sin 0 O cos 0 O 


—d p cos a 0 

(a,,, — d p sin a 0 ) cos 0 O 
-(am - dp sin a 0 ) sin 0 O 


0 0 0 


1 


(9) 


Here a 0 and 0 O are the normal pressure angle and the helix angle of the rack-cutter; a m 
is the half of the tooth width of the rack-cutter on middle line m — m (fig. 3), where 


7T 

4R 


( 10 ) 
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and P n is the normal diametral pitch of the rack-cutter, d p is the distance between middle 
line m-mof the rack-cutter and the origin Ob along axis x b , fig. 5. Parameter d p can be 
controlled to adjust the location of the contact path on the gear tooth surface. 

Surface S< of the gear rack-cutter is a plane that is represented in Sb as 

r?» = [u, 0 l,] T (") 


where (u t , l t ) are the surface parameters. 

Rack-cutter surface E t is represented in coordinate system S t by the matrix equation 


r,(u = M, t rf 


( 12 ) 


Equations (9), (11) and (12) yield 


r t 


(u t - d p ) COS Qo 

[(u t - d p ) sin cl 0 + a m ] cos 0 O + l t sin / 3 0 
— [(u t — d p ) sin a 0 + a m ] sin fi 0 + It cos 


(13) 


The unit normal to E t is represented in S t by equations 


_ N ‘ 

“ IN, | ’ 



dr t 

du t 


( 14 ) 


that yield 
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(15) 


n, 


[- 


sin a 0 cos a 0 cos 0 O — cos oc 0 sin j3 0 


Pinion Rack-Cutter Surface E c 

Rack-cutter E c generates the pinion. The normal section of rack-cutter surface E c (fig. 5) 
is a parabolic curve. We remind that the normal section of rack-cutter surface E t is a straight 
line directed along axis x b in fig. 5. The parabolic curve is in tangency with the zj-axis at 
point N(O b ). Rack-cutter surfaces E c and E t are in tangency along a straight line that is 
parallel to axes z a and z b and passes through point O b that coincides with point N . The 
deviation of the parabolic curve from the Zf,-axis affects the dimensions of the instantaneous 
contact ellipse. 

Rack-cutter surface E c is represented in S b as follows 



U c 


-a c Uc 



(16) 


where a c is coefficient of the parabolic normal section, and (u c , l c ) are the surface parameters 
of E c . 

Rack-cutter surface E c is represented in coordinate system S c by the matrix equation 

r c (u c ,/ c ) = M e6 r^ (17) 


Equations (9), (16) and (17) yield 
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(u c - d p ) cos a 0 + a c ul sin a Q 


[(u c — d p ) sin a 0 + a m ] cos 0 O — a c ul cos a 0 cos 3 0 + l c sin 0 O 


— [(u c — dp) sin a 0 + a m ] sin (d 0 + a c u 2 c cos a 0 sin 0 O + h cos p 0 J 


(18) 


The unit normal of E c is represented as 


n c 



w: 

du 


dv c 

W c 


7 % 

die' 


Equations (18) and (19) yield 


1 

e “ (1 + 4a2 u 2)0-5 


sin a 0 — 2 a c u c cos a„ 
—(cos a 0 + 2a c u c sin or 0 ) cos /3 0 
(cos a 0 + 2a c u c sin a„) sin /3 0 


(19) 


( 20 ) 


Using equations (13), (15), (18) and (20), it is easy to verify that surfaces E c and E t are 
in tangency along the axis when u c = u t = 0. 

4. Pinion and Gear Surfaces Generated by Rack-Cutters 

In the process for generation the two rigidly connected rack-cutters perform translational 
motion while the pinion and the gear perform rotational motions as shown in fig. 6. To 
provide a predesigned parabolic function of transmission errors for each cycle of meshing, it 
is necessary to observe certain relations between the motions of the rack-cutters and gears, 
respectively. 
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The angle V>pe of pinion rotation and the displacement s c of rack-cutter E e are related by 
the following linear function 


i>vc = - (2D 

T v 

Here: r p is the radius of the pinion pitch cylinder. 

The angle V> Jf of gear rotation and the displacement s t of rack-cutter E t are related as 
follows 


= ( 22 ) 

Here: N p and N g are the tooth numbers of the pinion and gear, respectively, and 0^ is 
the initial position angle of the gear for the modification gear rotation. 

Equation of Meshing between Rack-Cutter E c and Pinion E p 

The equation of meshing between rack- cutter E c and the pinion tooth surface E p is 
represented as 


/(«„ = N<=> • = 0 


( 23 ) 


where xj) pc is the angle of rotation of the pinion in the process for generation. The normal 
N(<=) to E c in S c can be obtained by equation (20), and the relative velocity of the pinion 
with respect to E c may be represented as 
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( 24 ) 


v (pc) _ w (p) x + T c ) - ( 0 T-p 0 ) 


1 

dt 


pc 


Here: R, = (UJT c )c = (r p r p ij) pc 0) r , r p = £ pc is the radius of the pitch cylinder of 
the pinion (fig. 6), = w^( 0 0 1 ) T . 

Substitution of equations (19) and (24) into (23), yields the following equation of meshing 
between E p and E c 


/(u c , / c , t/> pc ) = l e sin P 0 + r p xl) pc + a m cos 0 O + 


[(u e - dp) + 2 g^Ug] cos 0 O 
sin a„ — 2a c tt c cos a 0 


= 0 


(25) 


Surface of Pinion E p 

In the process of generation of pinion surface, rack-cutter E e perforins uniform translation 
and the workpiece performs uniform rotation (fig. 6(b)). The transformation matrix from 
system S c to S p can be represented as 


M pc = 


cost /> pc sin t/y 0 r p cos ^ + r p ip pc sin t/^ 

- sin xl> pc cos xp pc 0 -r p sin t/> pc + r p t/> pc cos ij) pc 

0 0 1 0 

0 0 0 1 


Pinion surface E p in system 5 P is represented as 


(26) 


r p (ti e ,/ c ,t/> pc ) = M pc r c 


(27) 


= -{-T^ — ^p) + 2a c u e , cos ^ + rp ^ pc + flm cos y3 0 }/ s in 
sma 0 — 2a c u c cosa 0 
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Substituting equation (26) into (27), we obtain equation of E p as 

r P (u c , 0 pc ) = T p (u p , dp) (28) 

Equation of Meshing between Rack-Cutter and Gear 

The equation of meshing between rack-cutter E t and the gear tooth surface E g is repre- 
sented as 

/(«., V,.) ■= n<‘> • v< a 'l = 0 (29) 

where rp gt is the angle of rotation of the gear in the process for generation. The unit normal 
n<^ to E t in S t is represented by equation (15), and v\ 9t ^ is the relative velocity of the gear 
with respect to rack-cutter E t . 

We recall that the rack-cutter E t performs translation with constant velocity, but the 
gear performs rotation with variable angular velocity that is represented as (see equation 
( 22 )) 

_ j 

«?>- o 0 - - 2a«.„c - 4.<°>)]j (30) 

The relative velocity vj 3 ^ is represented as 

v<*° = x(B, + r,)-(0 ^r„ 0) T ^p (31) 
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where 


R, = ( -r, 0 f (32) 

and T g = E gt is the radius of pitch cylinder of the gear (fig. 6). 

Substitution of equations (15), (30), (31) and (32) into (29) yields the following equation 

of meshing between E 3 and E t 


N, 


f(u t , It, $ g t) — ( u t — dp) cos (} 0 + sin a 0 (/t sin / 3„ + -jj-rgippc + flm cos 0 O ) 


2 aN g r g (if) pc - V> (0) ) a 

+ ; . tt^t— cos a 0 cos p 0 

Np — 2aNg(rp pc — V»<°>) 


(33) 


where 


</>», = - <■(-/>„= - 


(34) 


Surface of Gear E 3 

It must be remembered that the gear with the tooth surface E 3 performs rotation about 
its axis with varied angular velocity while rack-cutter E t performs uniform translation (fig. 
6(a)). The transformation matrix from system S t to S g can be represented as 
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M gt = 


— COS tfg gt 


— sin ij> g t 


0 

0 


sin i l> gt 0 

— cos ip g t 0 

0 1 

0 0 


N 

r g cos i/) gt + rgjf-rj>pc sin V> s t 

. . k, , 

Tg sin ipgt - rg—ippc cos ip gt 

lVg 

0 

1 


(35) 


Gear tooth surface E 3 in system S g can be represented as 


xj) gt ) — 


k = -{(u t - + (jrrg'l’pc + °m cos &)-r * 

sm a B I\g sin p 0 


1 


(36) 


2 aNgTjf^ ~ jff) 
AT p -2aiV 3 (^ c -0( 0 )) 


cot a 0 cot /? 0 } 


The derivation of equation (36) is based on transformation of equations (33) and (34). 
Equations (35) and (36) enable to represent the gear tooth surface in two-parameter form 
as follows 


r g{ u ti' t Pgt) — r g{ u gi@g) 


(37) 


5. Computerized Simulation of Meshing and Contact of Pinion-Gear Tooth 
Surfaces 

We consider that the surfaces of the pinion and the gear generated by worms Y, w and 
are represented in coordinate systems S p and S g , respectively. The fixed coordinate system Sf 
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is rigidly connected to the housing of the gear drive (figs. 7 and 8). The movable coordinate 
systems S p and S g are rigidly connected to the pinion and the gear, respectively. An auxiliary 
coordinate system is applied for simulation of meshing when the gear axis is crossed or 
intersected with the pinion axis instead of being parallel, and when the shortest distance 
between the pinion and gear axes is changed. The errors of misalignment are referred to the 
gear. The misalignment angle A7 is decomposed into two components, A 7* and A~f v that 
represent the crossing angle and the intersection angle, respectively. The pinion performs 
rotational motion about the z/-axis. The axis of gear rotation is z^. The shortest distance 
between the axes of rotation is designated as Epg. 

The rotation matrices from system to Sf for crossed and intersecting angles are rep- 
resented in the followings (fig. 8) 


L f h 


-10 0 
0 — cos A7 x sin A y x 

0 sin A7 x cos A-y x 


( 38 ) 


L fh = 


— cos A7 v 0 — sin A7 V 

0-10 

— sin A7 v 0 cos A7 V 


( 39 ) 


We represent the pinion and gear tooth surfaces, E p and E 3 , and their unit normals in 
coordinate system Sf. The conditions of continuous tangency of surfaces E p and E 3 are 
represented by the following equations [1,2]. 


'/’Vp^p^p) = r /' ) KA>^) 


( 40 ) 
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n / P) ( U PT 9pi <Pp) = ^{Ug, dg, <j)g ) 


( 41 ) 


Vector equation (41) provides only two independent equations since |njr | = |n^| = 1. 
The total number of independent equations provided by (40) and (41) is five that relate six 
parameters 


/«( U pj dp, <f>p , Ug , 6g , <f>g) 0 


(i = 1,2, ...5) 


(42) 


The continuous solution of the system of nonlinear equations (42) is based on the following 
procedure: 

(1) Using an initial guess, we determine a set of parameters that satisfy equation system 
(42). Thus 


p {0) = (4 MVM°US ,) ) 


(43) 


(2) One of the variable parameters, say <f> p , is chosen as the input one, and is supposed 
that the Jacobian 


D{Up, 6p , Ug , dg , <t>g) 


differs from zero. The derivatives in the Jacobian are taken at point P^°\ 

(3) Then, equation system (42) can be solved in the neighborhood of P ^ by functions 
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<t>g{<t>p), U p(<M’ ^p(^p)> U 5 (<Ap)’ 


(45) 


(4) Vector function r p (u p ,0 p ) that determines the pinion surface E p and functions u p (<£ p ), 
^p(^p) enable to determine the path of contact on E p . 

(5) Similarly, we can obtain the path of contact on the gear surface E 3 using vector 
function r g {u g ,6 g ) and functions u g (4> ? ), 9 g {4> g ). 

(6) The paths of contact on pinion and gear tooth surfaces slightly deviate from helices 
in the case of an aligned gear drive. The line of action for an aligned gear drive (the set of 
contact points in Sf) slightly deviates from a straight line that is parallel to the gear axes. 

(7) The transmission function <f> g (<t> p ) deviates from the ideal transmission function, and 
the function of transmission errors coincides with the predesigned parabolic function. 

(8) The determination of dimensions and orientation of the instantaneous contact ellipse 
needs the knowledge of the principal curvatures and directions of contacting surfaces and the 
elastic approach of surfaces. This problem can be substantially simplified if the pinion-gear 
principal curvatures and directions are expressed in terms of the principal curvatures and 
directions of the generating surfaces and parameters of motion [1 ,2]. 

6. Numerical Example 

The method developed in this report is illustrated with the example discussed below. The 
design parameters of the pinion and gear are listed in Table 1. The numerical simulation 
of meshing is performed for an aligned and misaligned gear drives with various errors of 
alignment for the pinion and gear. 

Case 1. Aligned gear drive 

Figure 9 shows the transmission errors for the aligned gear drive. The TCA performed 
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Table 1: Design parameters of pinion and gear 



pinion 

gear 

tooth number 

N p = 20 

N g = 100 

normal diametral pitch 

Fn = 5 - 

Pn = 5 — 


in 

in 

normal pressure angle 

P 

0 

II 

to 

0 

0 

a 0 = 20° 

helix angle on pitch cylinder 

0 

0 

CO 

II 

II 

CO 

0 

0 

tooth length 

L=1.6 in. 

L=1.6 in. 

modification coefficient 

a c = 0.0008 

a - 0.0014 

elastic approach 

& = 0.007 mm. 

8 = 0.007 r 


confirms that the predesigned parabolic function of transmission errors exists. The maxi mum 
transmission error is approximately 8 arc seconds. Figure 10 shows the contact pattern and 
the contact path. The path of contact on the tooth surface is in the longitudinal direction. 
The major axis of the instantaneous contact ellipse is 6 mm for the assumed elastic approach 
of the surfaces equal to 0.007 mm. 

Case 2. The pinion-gear rotation axes are crossed 

Figures 11 and 12 show the transmission errors and the contact pattern for the case when 
the crossing angle Across is 4 arc minutes. The maximum transmission error is 8 arc seconds 
and contact paths are shifted up and down on the gear and pinion surfaces, respectively. 
Figures 13 and 14 show the transmission errors and contact pattern for A7 CT0SS = —4 arc 
minutes. 

Case 3. The pinion-gear rotation axes are intersected 

Figures 15 and 16 show the transmission errors for the misalignment Winter sect = 4 
arc minutes. Figures 17 and 18 show the transmission errors and contact pattern for the 
mentioned above error of alignment. 
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Case 4. Influence of error of lead angle 

Figures 19 and 20 show the transmission errors and contact pattern when the error of 

the lead angle is 4 arc minutes. 

Figures 21 and 22 show the transmission errors and contact pattern when the error of 
the lead angle is —4 arc minutes. 

For all of above cases, the maximum transmission error does not exceed 8 arc second 
(with very small deviations of this value). 

7. Conclusion 

From the analytical study presented in this report the following conclusions can be drawn: 

(1) The interaction of a parabolic and a linear functions of transmission errors has been 
discussed to prove the possibility to absorb almost, linear functions of transmission 
errors caused by misalignment. 

(2) Mismatched surfaces of two rack-cutters for generation of modified involute gears have 
been proposed. 

(3) Generation and geometry of pinion-gear modified tooth surfaces have been determined. 

(4) Computerized simulation of meshing and contact of pinion-gear tooth surfaces has been 
developed. 

(5) An algorithm for determination of relations between the curvatures of the generating 
and the generated surfaces has been developed. 

(6) An algorithm for determination of the contact ellipse has been developed. 

(7) Directions for users of application of developed computer programs for the design of 
gears with the modified geometry and computerized simulation of their meshing and 
contact have been developed (Appendix C). 
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Fig. 3 Normal sections of rack-cutters 



25 


Parallel to gear axes 



Fig. 4 Orientation of rack-cutters with respect to gear axes 
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(a) 




Fig. 8 Coordinate systems and Sf 
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TRANSMISSION ERRORS (second! 


(NO ASSEMBLY ERRORS! 



PINION ROTATION ANGLE (degree) 


Fig. 9 Transmission errors for the aligned gear drive 
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Fig. 10 Contact paths and pattern for the aligned gear drive 
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CROSSING ANGLE BETWEEN RXES: 4 MINUTES 



Fig. II Transmission errors for the alignment error 
A7 cross = 4 arc minutes 
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CROSSING RNGLE BETWEEN RXES: 4 MINUTES 
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Fig. 12 Contact pattern for the alignment 
A7 cr0 ss = 4 arc minutes 


error 
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TRANSMISSION ERRORS I second! 


CROSSING ANGLE BETWEEN AXES: -* MINUTES 



Fig. 13 Transmission errors for the alignment error 
Across = - 4 arc minutes 


35 


CROSSING ANGLE BETWEEN AXES: MINUTES 



Fig. 14 Contact pattern for the alignment error 
A 7 cr o Sg = - 4 arc minutes 
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INTERSECTING ANGLE EETWEEN AXES: 4 MINUTES 



PINION ROTATION ANGLE (degree) 


Fig. 15 Transmission errors for misaligned gear drive with 
A7 intersectting = 4 arc minutes 


37 



INTERSECTING RNGLE BETWEEN RXES: 4 MINUTES 



AXIS Z MM 


Fig. 16 Contact pattern and paths for the misaligned 
gear drive with A7 iatersecttin5 = 4 arc minutes 
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INTERSECTING ANGLE BETWEEN AXES: -+ MINUTES 



Fig. 17 Transmission errors for the misaligned 
gear drive with A7 intersectUng = - 4 arc minutes 
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INTERSECTING HNGLE BETWEEN RXES: MINUTES 
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Fig. 18 Contact pattern for the misaligned 
gear drive with A7 i[ller « cttiri! = - 4 arc minutes 


40 



ERROR OF LERD ANGLE: 4 MINUTES 



PINION ROTRTION ANGLE (degree) 


Fig. 19 Transmission errors for the misaligned 
gear drive with \\ 0 = 4 arc minutes 
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ERROR OF LERD RNGLE: 4 MINUTES 



AXIS Z MM 


Fig. 20 Contact pattern and contact path for misaligned 
gear drive with AX 0 = 4 arc minutes 
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PINION ROTRTION ANGLE (degree) 


Fig. 21 Transmission errors for the misaligned 
gear drive with AXq = ^ arc m i nu tes 
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Appendix A. Relations between the Curvatures of the Generating and Generated 
Surfaces 

Direct Relations between Principal Curvatures and Directions of Mating Sur- 
faces 

The main advantage of this approach (proposed by Litvin) is the possibility to determine 
the principal curvatures and directions of the generated surface in terms of principal curva- 
tures and directions of the generating tool surface, and the parameters of motion. In this 
case, the tool surfaces for the generation of the gear and the pinion tooth surfaces are rack 
cutters. The equations developed permit a simplified computational procedure. 

The system of equations that relate the principal curvatures and directions of the gen- 
erating and generated surfaces is as follows. Consider that unit vectors e f and e* represent 
the principal directions on the tool surface Ei at point P of tangency of surfaces E t and E 2 
(fig. Al). The principal curvatures on the mating surfaces K f and K h of the tool are given; 
the parameters of motion (see below) are also given. 

The goal is to determine angle a that is formed by unit vectors e/ and e 3 , and principal 
curvatures k s and K q . (The unit vectors e, and e, represent the principal directions on 
surface E 2 ). The system of equations for determination of <r, k s , and K q is as follows. 


tan 2 a = 


2bi5& 2 5 

b 2S K h)*33 


(Al) 


2 & 15 6 2 5 *> 2 5 ~ ^15 ~ (*/ * 0^33 

f 33 sin 2 a *33 cos 2 <r 


(A2) 


K q + K, = Kf + K h + 


^15 ^25 

*33 


Here: 


(A3) 
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bis = -M 12 > • e fc ) - Mv ( 12) ' e/) 


m 


6 25 = M 12) • e/) - K fc (v (12) ■ e*) 


(AS) 


<33 = -n • [(w (1) x v t (2) ) - (u>< 2 ) x vj?)] 

+(n x u>( 12 >) • v( 12 ) — «/(v( 12 ) • e/) 2 — k*(v( 12 ) • e,,) 2 
The nomenclature for equations (A4) to (A6) is described as follows: 

u/ 1 ) angular velocity of the generating tool 

0 ,( 2 ) angular velocity of the generated gear 

o;( 12 ) defined as u/ 1 ) - o;( 2 ) 

v[*) transfer motion velocity of the generating tool 

vj 2 ^ transfer motion velocity of the generated gear 

v (12) defined as v]^ - v t (2) 

n surface unit normal vector 


(A6) 


The equations discussed above are used in the TCA program for determination of the 
contact ellipse at the points of contact path of the modified helical gear drive. 


Numerical Example 

The input and output for the determination of the principal curvatures of the pinion 
tooth surface are shown in Tables A1 and A2. The input and output for the determination 
of the principal curvatures of the gear tooth surface are shown in Tables A3 and A4. 
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Table A1 Input Data 


Description 

Symbol 

Values 

first principal curvature of tool 

K, 

-0.0016 (1/mm) 

second principal curvature of tool 

mm 

0.0 (1/mm) 

first principal direction of tool 

BUM 

[0.0016 - 0.4999 0.8660] r 

second principal direction of tool 

e* 

[0.9387 0.2992 0.1710] T 

angular velocity of tool (1/sec) 


[0.0 0.0 0.0] r 

angular velocity of pinion (1/sec) 


M'lflf'IM 

transfer velocity of tool 
(mm/ sec) 

vS> 

[0.0 58.6588 0.0] r 

transfer velocity of pinion 
(mm/ sec) 

vl 2) 

[-0.1514 58.5951 0.0] T 

surface normal of tangent point 

n 

[-0.3420 0.8138 0.4698] T 


Table A2 Output Data 


Description (for pinion) 

Symbol 

Values 

first principal curvature 

K,, 

-0.001543 (1/mm) 

second principal curvature 

K„ 


first principal direction 

e. 

[0.1751 - 0.4361 0.8827] T 

second principal direction 

e ? 

[0.9222 0.3865 0.0080] 7 
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Table A3 Input Data 


Description 

Symbol 

Values 

first principal curvature of tool 

Kf 

0.0 (1/mm) 

second principal curvature of tool 


0.0 (1/mm) 

first principal direction of tool 

e / 

[0.9397 0.2962 0.1710] r 

second principal direction of tool 

e* 

[0.0 - 0.5 0.8660] t 

angular velocity of tool (1/sec) 


[0.0 0.0 0.0] T 

angular velocity of gear (1/sec) 



transfer velocity of tool 
(mm/sec) 

vg> 

[0.0 58.6588 0.0] T 

transfer velocity of gear 
(mm/sec) 

vi? 

[-0.2096 58.7218 0.0] T 

surface normal of tangent point 

n 



Table A4 Output Data 


Description (for gear) 

Symbol 

Values 

first principal curvature 



second principal curvature 

*7 

0.0 (1/mm) 

first principal direction 

e. 

[0.9219 0.3874 0.0] r 

second principal direction 

e ? 

[0.1820 - 0.4331 0.8827] T 
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Fig. A1 Principal Directions 
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Appendix B. Contact Ellipse 


Determination of Dimensions and Orientation of Instantaneous Contact Ellipse 

The gear tooth surfaces are in point contact at every instant. Due to elastic deformation 
of gear tooth surfaces the contact is spread over an elliptical area and the center of the ellipse 
coincides with the instantaneous contact point. The hearing contact is formed as the set of 
instantaneous contact ellipses. 

The dimensions and orientation of the instantaneous contact ellipse can be determined 
using the data about the principal curvatures and directions of the contacting surfaces, and 
the elastic approach of the surfaces. The elastic approach depends on the applied load but 
we will consider it as a given value that is known from experimental data. 

The determination of the instantaneous contact ellipse is based on the following equations 
(proposed by Litvin): 


where 


cos2a^ = 


sin2a^ = 


2a = 2 


9i — 92 cos 2 <t 

{ 9 i ~^ 9 i 92 cos2<r + ^) 1 /2 

g 2 sin 2<r 

(Si - 2 gi 92 cos 2 <t + gj ) 1/2 

1/2 


A 

1/2 

8 

\J 

, 26 = 2 


A 


~B 


A = ^[ 4 ^ - 4 2) ~ ( 9i ~ ^9192 cos 2 a + gl) 1/2 ] 
B = '-[k^ - 4 2) + 0 9 1 - 29x92 cos 2<r + g 2 ) 1/2 j 


(B 1 ) 
(£ 2 ) 
(B 3) 

(B 4) 
(B 5) 


„(>•) _ „(0 , -W _ _(p _ JO 

— K i + K n » 9% — K i K n 


(Be) 
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Here (fig. Bl) a (1) is the angle that is formed by axis rj of the contact ellipse with the unit 
vector e; ^ of the principal direction on surface Ej; cr is the angle formed by unit vectors 
and of the principal directions of the contacting surfaces; 2 a and 26 are the axes of 
the contact ellipse; 8 is the elastic approach; and and are two principal directions of 
tooth surface i. 


Numerical Example 

The input and output for the determination of the contact ellipse are shown in Tables 
Bl and B2. 


Table Bl Input Data 


Description 

Symbol 

Values 

pinion first principal curvature 

K* 

-0.001543 (1/mm) 

pinion second principal curvature 

K, 

0.03907 (1/mm) 

pinion first principal direction 


[0.1751 - 0.4361 0.8827] T 

pinion second principal direction 

e ? 

[0.9222 0.3865 0.0080] r 

gear first principal curvature 

K f 

-0.007836 (1/mm) 

gear second principal curvature 


0.0 (1/mm) 

gear first principal direction 

e / 

[-0.9219 - 0.3874 0.0] T 

gear second principal direction 

e* 

[-0.1820 0.4331 0.8827] T 

elastic approach 

6 

0.007 (mm) 


Table B2 Output Data 


Description 

Symbol 

Values 

long axis of contact ellipse 

2a 

6.026 (mm) 

short axis of contact ellipse 

26 

1.092 (mm) 

angle between long axis and 

the first principal direction of pinion 


89.87 (deg) 
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Fig. B1 Orientation and dimensions of contact ellipse 
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Appendix C. Directions for Users of Application of Computer Program 
C.l Introduction 

The name of the computer program is HELTCA.FOR. It is written in FORTRAN77 
language. The operating system is CMS-9.0. A subroutine DNEQNF to solve a system of 
nonlinear equations should be available in the Math-Library or working environment. The 
subroutine is not included in the program. The program will call the subroutine DNEQNF 
several times. 

C.2 Input Block 

The input block consists of three parts: (1) design parameters of pinion and gear; (2) the 
controlled modification parameters; and (3) parameters for TCA. 

Part 1 . Design parameters of pinion and gear 

In the beginning of the computer program, you can read the following lines: 

C All COEFFICIENT FOR TRANSFORMATION OF DEGREE TO RADIAN 

A1I=DAC0S(— l.D0)/180.D0 
C 

C... KHD=1 FOR RIGHT-HAND PINION AND LEFT-HAND GEAR 
C... KHD=2 FOR LEFT-HAND PINION AND RIGHT-HAND GEAR 

KHD=2 

If you write “KHD=1”, the computer will use the necessary equations for the case of 
right-hand pinion and left-hand gear. The computations will be accomplished for a left- 
hand pinion and right-hand gear if you use “KHD=2”. 

Then, the variable definition for the pinion follows: 

C 

C... INPUT THE DESIGN PARAMETERS OF PINION 
C... 

C... TNI GEAR NUMBER OF TEETH 

C... PN1 NORMAL DIAMETRAL PITCH (1/MM) 

C... PSIN1 NORMAL PRESSURE ANGLE (RAD.) 

C... BETAP1 HELICAL ANGLE ON PITCH CYLINDER (RAD.) 

C... ADG1 ADDENDUM (MM) 
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C... DEG1 DEDENDUM (MM) 

C... LAMDP1 LEAD ANGLE ON PITCH CYLINDER (RAD.) 

C... FW1 FACE WIDTH (MM) 

C... RPT1 RADIUS OF PITCH CYLINDER (MM) 

C... RBT1 RADIUS OF BASE CYLINDER (MM) 

C... RATI RADIUS OF ADDENDUM CYLINDER (MM) 

C... RDT1 RADIUS OF DEDENDUM CYLINDER (MM) 

C... PSIT1 TRANSVERSE PRESSURE ANGLE (RAD.) 

C... L AMDB1 ....LEAD ANGLE ON BASE CYLINDER (RAD.) 

C 

In accordance with our numerical example (see Table 1), the following data would be 


used: 


TN1=20.D0 

PN1=5.D0/25.4D0 

PSIN1=A11*20.D0 

LAMDP1=A11*60.D0 

BETAP1=A11*30.D0 

FW1=25.4D0*1.6D0 

ADG1=1.D0/PN1 

DEG1=1.25D0/PN1 


The following variables are used for the gear: 

C 

C... INPUT THE DESIGN PARAMETERS OF GEAR 
C 

C... TN2..............GEAR NUMBER OF TEETH 

C .. PN2 NORMAL DIAMETRAL PITCH (1/MM) 

C... PSIN2 NORMAL PRESSURE ANGLE (RAD.) 

C .. BETAP2 HELICAL ANGLE ON PITCH CYLINDER (RAD.) 

C... ADG2 ADDENDUM (MM) 

C... DEG2 DEDENDUM (MM) 

C... LAMDP2 LEAD ANGLE ON PITCH CYLINDER (RAD.) 

C... FW2 FACE WIDTH (MM) 

C .. RPT2 RADIUS OF PITCH CYLINDER (MM) 

C RBT2 RADIUS OF BASE CYLINDER (MM) 

C... RAT2 RADIUS OF ADDENDUM CYLINDER (MM) 

C RDT2 RADIUS OF DEDENDUM CYLINDER (MM) 

C PSIT2 TRANSVERSE PRESSURE ANGLE (RAD.) 

C... LAMDB2 LEAD ANGLE ON BASE CYLINDER (RAD.) 

C 


From the design parameters listed in Table 1 , we have 
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TN2=100.D0 

PN2=5.D0 /25.4D0 

All=DACOS(-l.D0)/180.D0 

PSIN2=A11*20.D0 

LAMDP2=A11*60.D0 

BETAP2=A11*30.D0 

FW2=25.4D0*1.6D0 

ADG2=1.D0/PN2 

DEG2=1.25D0/PN2 


If KHD=2, the computer program will change the values of some design parameters as 
follows 


IF(KHD.EQ.2) THEN 
LAMDP1=-LAMDP1 
BETAP1=-BETAP1 
LAMDP2= -L AMDP2 
BETAP2=-BETAP2 
ENDIF 

The computer program call the following subroutines “DATAT1” and “DATAT2” to 
calculate other tooth element proportions and output the whole data in file 55 (see below). 

CALL DATAT1 
CALL DATAT2 


Part 2. Control of modification parameter for application in the TCA program 
At this stage we can read 


C THE FOLLOWING DATA IS FOR THE TO-BE CONTROLLED MODIFICATION 

PARAMETERS 



AA... MODIFICATION PARAMETER OF GEAR 
AP... MODIFICATION PARAMETER OF PINION RACK-CUTTER 
DP TANGENT POINT N OF PROFILES OF PINION & GEAR RACK-CUTTERS 
THET2P. .INITIAL ANGLE FOR MODIFICATION OF GEAR (RAD.) 


You must input the four parameters for modification of pinion and gear surfaces, for 


example: 


AA=— 0.0014D0 
THET2P=-0.08D0 
AP=— 0.0008 

DP=— DSIN(PSIN1)*DCOS(-1.DO)*RPT1/TN1/8.DO 
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The above four controlling parameters should be tried several times is order to obtain 
better contact pattern and transmission errors optimal for a given design. 


Part 3. Parameters for TCA 

In this part, the alignment errors expected should be input: 

C 

C PARAMETERS FOR TCA 

C 

C.. KM .. SWITCH 1 FOR CROSSING ANGLE A 7ar & 2 FOR INTERSECTION 
ANGLE A 7„ 

C.. DGAM..ANGLE OF MISALIGNMENT(CROSSING OR INTERSECTION) (ARC 
MINUTE) 

C.. DEE... CHANGE OF CENTER DISTANCE (MM) 

C 

If a crossing angle of misalignment is considered, should input “KM=1”. Input “KM=2” 
if an intersection angle of misalignment is considered. For an aligned gear drive, input 
u DGAM=0.0”. For instance, if Aj r = 4', A E = 0, input the following lines: 

KM=1 

C KM=2 

DGAM=4.0* A1 1 /60.D0 
DEE=0.000D0 

Then you will read the following sentence: 

C 

C.... THE INPUT BLOCK IS READY 

i 

Usually, you cannot make changes anything after this step. 


C.3 Output Block 

After the input block is filled out, you can read the following explanation for the output 


files: 

C 

C... OUTPUT DATA FILES ARE THE FOLLOWINGS 
C 

C... FILE 55... TOOTH PROPORTIONS OF PINION AND GEAR 
C... FILE 85... CONTACT PATH ON PINION SURFACE (2D) 

C... FILE 86... CONTACT PATH ON GEAR SURFACE (2D) 

C... FILE 87... LENGTH AND DIRECTIONS OF CONTACT ELLIPSE ON PINION 
AND GEAR SURFACES (2D) 
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C... FILE 90... TRANSMISSION ERRORS 
C 

C 

C 

1) File 55 

In File 55 the information about the pinion and gear are listed. 

2) Files 85 and 86 

There are two coordinates in File 85 for each contact point of the pinion: Radial (£p+y p ) 0 ' 5 
and axial z p . There are two coordinates in File 86 for each contact point of the gear: Radial 
{x] + y*) 0 5 and axial z g . 

3) File 87 

There are 5 values in File 87 for each pair of contact points of the pinion and gear. The 
first one is the value of the major semi- axis of the contact ellipse. The second and third 
values are the cosine directions of the major axis of the contact ellipse on the pinion tooth 
surface. The last two values are the cosine direction of the major axis of the contact ellipse 
on the gear tooth surface. 

C.4 Computer Program 
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c 


TCA FOR MODIFIED HELICAL GEARS 

PROGRAM HELTCA 
IMPLICIT REAL*8(A-H,0-Z) 

REAL*8 XI (9) , X(9) , F(9) 

REAL* 8 LF1 (3 , 3) ,LH2(3,3) ,L1F(3,3) ,L2H(3,3) 

REAL* 8 RIF (3) , R2F ( 3 ) ,N1F(3) , N2F ( 3 ) 

REAL* 8 DPHI (2 , 180) , DDPHI (2,180,4) 

REAL* 8 ELAL(180) , ELI (2, 180) ,EL2(2,180) 

REAL* 8 LFH (3,3), LHF (3,3) 

REAL* 8 VI (3) , V2 ( 3 ) , V3 ( 3 ) ,V4(3) ,V5(3) 

REAL* 8 RG2 ( 3 ) , NG2 ( 3 ) 

REAL*8 RG1(3) , NG1(3) 

REAL*8 UI ( 3 ) , UJ ( 3 ) , UK ( 3 ) 

REAL* 8 LAMDP2 , LAMDB2 
REAL* 8 LAMDP1 , LAMDB1 
REAL*8 UT2 ,KT21,KT22 

REAL*8 AVC1 ( 3 ) , VTR1 ( 3 ) , AVC2 ( 3 ) , VTR2 ( 3 ) 

REAL*8 KSIG1 , KSIG2 , KFF , KHH 

REAL* 8 EFN(3) ,EHN(3) ,W1VT2(3) ,WV12(3) ,W2VT1(3) , KF , KH , KS , KQ 
REAL* 8 KM2 ,KT2 
REAL*8 KFP , KHP , KFG , KHG 
COMMON /A300/ ES(3) ,EQ(3) 

COMMON /A310/ KFF, KHH 
COMMON /A340/ EFF(3) ,EHH(3) 

COMMON /A360/ A,B, SI (3) , FI (3 ) 

COMMON/ A200/ Wl(3) ,W2(3) ,W12(3) ,VT1(3) ,VT2(3) ,V12(3) 

COMMON /A210/ EX(3) ,EF(3) ,EH(3) 

COMMON /A212/ EF2 (3 ) , EH2 (3) , KF2 , KH2 ehe i-sv 

COMMON /A220/ KF,KH 
COMMON /A230/ ET(3),EM(3) 

COMMON /A380/ KS,KQ 

COMMON /A400/ VT11(3) ,VT12 (3) ,VT21(3) ,VT22 (3) 

COMMON /A401/ KHP , KFP , KHG , KFG 
EXTERNAL FCNG , FCNC , FCNT 
COMMON /AXIS/ UI , UJ , UK 
COMMON /NET/ RR,DD 
COMMON /DD/ DF , KPRI 

COMMON /DATT2 / TN2 , PN2 , PSIN2 , BETAP2 , ADG2 , DEG2 , LAMDP2 , 

& UP2 , FW2 , RPT2 , RBT2 , RAT2 , RDT2 , PSIT2 , LAMDB2 
COMMON /DATT/ TN1,PN1,PSIN1,BETAP1,ADG1,DEG1,LAMDP1, 

& FW1 , RPT1 , RBT1 , RATI , RDT1 , PSIT1 , LAMDB1 
COMMON /B2 /XNP1 , YNP1 , ZNP1 , XNP2 , YNP2 , ZNP2 
COMMON /B4 /XI , Y1 , Z1 , XN1 , YN1 , ZN1 
COMMON /B5/X2,Y2,Z2,XN2,YN2,ZN2 
COMMON /B6 /THET2P , DGPHI2 

COMMON /Wl/ ETAW1 , UPP1 , SPP1 , ETAW2 , UPP2 , SPP2 
COMMON /SGI/ RG1 , NG1 , AP , DP 
COMMON /SG2/ RG2,NG2,AA 

COMMON /ATT/ PHI1,PHI2,R1F,R2F,N1F,N2F,LFH,CC,DGAM,DPHI2 
COMMON /MVT/ LF1 , LH2 
COMMON /AST/ ICONT 
COMMON /ATS/ DPHI1 
C... All. ... .COEFFICIENT FOR TRANSFORMATION OF DEGREE TO RADIAN 
All=DACOS (-1 . DO ) / 180 . DO 

C... KHD=1 FOR RIGHT-HAND PINION AND LEFT-HAND GEAR 
C. . . KHD=2 FOR LEFT-HAND PINION AND RIGHT-HAND GEAR 
KHD=2 

C... INPUT THE DESIGN PARAMETERS OF PINION 


C. .. TNI GEAR NUMBER OF TEETH 

C... PN1 NORMAL DIAMETRAL PITCH (MM) 
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C... PSIN1 NORMAL PRESSURE ANGLE (RAD.) 

C... BETAP1 LEADING ANGLE OF THE HELIX ON PITCH CYLINDER (RAD.) 

C... ADG1 ADDENDUM (MM) 

C. . . DEG1 DEDENDUM (MM) 

C... LAMDP1. .. .HELIX ANGLE ON PITCH CYLINDER (RAD.) 

C... FW1 FACE WIDTH (MM) 

C... RPT1 RADIUS OF PITCH CYLINDER (MM) 

C... RBT1 RADIUS OF BASE CYLINDER (MM) 

C... RATI RADIUS OF ADDENDUM CYLINDER (MM) 

C. .. RDT1 RADIUS OF DEDENDUM CYLINDER (MM) 

C... PSIT1 TRANSVERSE PRESSURE ANGLE (RAD.) 

C... LAMDB1 HELIX ANGLE ON BASE CYLINDER (RAD.) 


TN1=20.D0 
PN1=5.D0/25.4D0 
PSIN1=A11*20 .DO 
LAMDP1=A11*60 . DO 
BETAP1=A11*3 0 . DO 
FW1=25 . 4D0*1 . 6D0 
ADG1=1 . DO/PNl 
DEG1=1 . 25D0/PN1 

C. . . INPUT THE DESIGN PARAMETERS OF GEAR 


C... TN2 GEAR NUMBER OF TEETH 

C... PN2 NORMAL DIAMETRAL PITCH (1/MM) 

C. .. PSIN2 NORMAL PRESSURE ANGLE (RAD.) 

C. .. BETAP2 LEADING ANGLE OF THE HELIX ON PITCH CYLINDER (RAD.) 

C... ADG2 ADDENDUM (MM) 

C... DEG2 DEDENDUM (MM) 

C. .. LAMDP2 HELIX ANGLE ON PITCH CYLINDER (RAD.) 

C. . . FW2. FACE WIDTH (MM) 

C... RPT2 RADIUS OF PITCH CYLINDER (MM) 

C. .. RBT2 RADIUS OF BASE CYLINDER (MM) 

C. .. RAT2 RADIUS OF ADDENDUM CYLINDER (MM) 

C... RDT2 RADIUS OF DEDENDUM CYLINDER (MM) 

C... PSIT2 TRANSVERSE PRESSURE ANGLE (RAD.) 

C... LAMDB2 HELIX ANGLE ON BASE CYLINDER (RAD.) 


TN2=100 . DO 
PN2=5.D0/25.4D0 
All=DACOS (-1.D0) /180.D0 
PSIN2=A11*20 . DO 
LAMDP2=A11*60.D0 
BETAP2=A11*30.D0 
FW2=25 . 4D0*1 . 6D0 
ADG2=1 . DO /PN2 
DEG2=1 . 25D0/PN2 
IF(KHD.EQ.2) THEN 
LAMDP1=-LAMDP1 
BETAP1=-BETAP 1 
LAMDP2=— LAMDP2 
BETAP2=-BETAP2 
ENDIF 

CALL DAT ATI 
CALL DATAT2 

C... THE FOLLOWING IS FOR CONTROLLING PARAMETERS 

C... AA. .. MODIFICATION PARAMETER OF GEAR 

C... AP... MODIFICATION PARAMETER OF PINION RACK-CUTTER 

C. . DP TANGENT POINT N OF PROFILES OF PINION £ GEAR RACK-CUTTERS 

C. . THET2P. .INITIAL ANGLE FOR MODIFICATION OF GEAR (RAD.) 

C 

C PARAMETERS FOR TCA 
C 
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c. . 
c. . 
c. . 


KM . . SWITCH 1 FOR CROSSING & 2 FOR INTERSECTING MISALIGNMENT 
DGAM. .ANGLE OF DISALIGNMENT (CSOSSING OR INTERSECTING) (ARC MINUTE) 
DEE. . . CHANGE OF CENTER DISTANCE (MM) 

KM=1 
C KM=2 

DGAM=0 . DO 
DEE=0. 000D0 
AA=-0. 00 14 DO 
THET2P=-0 . 08D0 
AP=-0. 0008 

DP=-DSIN (PSIN1) *DCOS (-1 . DO) *RPTl/TNl/8 . DO 
C... THE INPUT BLOCK IS OVER HERE 
C 

C... OUPTPUT DATA FILES ARE THE FOLLOWINGS 
C 

C... FILE 55... TOOTH PROPORTIONS OF PINION AND GEAR 
C... FILE 85... CONTACT PATH ON PINION SURFACE (2D) 

C... FILE 86... CONTACT PATH ON GEAR SURFACE (2D) 

C... FILE 87... DIRECTIONS OF LONG AXIS OF CONTACT ELLIPS (2D) 

C... FILE 90 ...TRANSMISSION ERRORS 
DO 901 1=1,3 
UI (I) =0 .DO 
UJ (I) =0 .DO 
UK(I) =0 . DO 
901 CONTINUE 
UI (1) =1 • DO 
UJ (2) =1 . DO 
UK(3 ) =1 . DO 

C .. EE2. . . . . GEAR RATIO 
EE2=TN1/TN2 

C .. CC CENTER DISTANCE OF GEAR DRIVE 

CC=RPT1+RPT2+DEE 

C . . CALCULATE CONTACT POINT ON MEAN SECTION WITHOUT MISALIGNMENT 
ICONT=l 
N=6 

ERRREL=0 . ID-6 
ITMAX=1000 
C 

CALL MIAL(KM, 0.D0,LFH) 

C 

XI(1)=0.D0 
XI (2) =0 .DO 
XI(3)=0.D0 
XI(4)=0.D0 
XI (5) =0. DO 
XI (6) =0. DO 
DD=0 . DO 
C 

CALL DNEQNF ( FCNG , ERRREL , N , ITMAX , XI , X , FNORM) 

C 

PHISS1=X(6) 

PHI1MEA=X ( 6 ) 

C .. CALCULATE CONTACT POINT ON EDGE SECTION WITHOUT MISALIGNMENT 
DD=-0.5D0*FW1 
ERRREL=0 . ID-8 
c 

CALL DNEQNF (FCNG, ERRREL, N, ITMAX, XI, X, FNORM) 

C 

PHISS2=X (6) 

C .. CALCULATE CONTACT POINT ON EDGE SECTION WITH MISALIGNMENT 
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N=6 

ERRREL=0 . ID— 6 
ITMAX=1000 

C 

CALL MIAL(KM,DGAM,LFH) 

C 

XI(1)=X(1) 

XI ( 2 ) =X ( 2 ) 

XI ( 3 ) =X ( 3 ) 

XI ( 4 ) =X ( 4 ) 

XI ( 5 ) =X ( 5 ) 

XI (6)=X(6) 

C 

CALL DNEQNF ( FCNG , ERRREL , N , ITMAX , XI , X , FNORM) 

C 

1330 CONTINUE 
PHI1=X(6) 

PHI 1STA=X ( 6 ) 

C . . THE FOLLOWING IS FOR TCA 
ICONT=2 
N=5 

ERRREL=0 . ID-5 
ITMAX=400 

STP= (PHISS2-PHISS1) /36 . DO 
NN=72 

DO 1010 1=1, NN 

PHI1=PHI1STA- (1-1) *STP 
XI(1)=X(1) 

XI ( 2 ) =X ( 2 ) 

XI(3)=X(3) 

XI(4)=X(4) 

XI ( 5 ) =X ( 5 ) 

C 

CALL DNEQNF (FCNG, ERRREL, N, ITMAX, XI, X, FNORM) 

IF(RG1(3) .GT. (0.5D0*FW1) ) GOTO 1011 

C... ROTATING VELOCITY OF CUTTER IN CUTTER SYSTEM 
W1(1)=0.D0 
W1(2)=0.D0 
W1 (3 ) =0 .DO 

C... ROTATING VELOCITY OF PINION IN CUTTER SYSTEM 
W2 (1) =0 . DO 
W2(2)=0.D0 
W2(3)=1.D0 

C... NORMAL OF CUUTER IN CUTTER SYSTEM 
EX(1)=XNP1 
EX ( 2 ) =YNP1 
EX ( 3 ) =ZNP1 

C. . . TRANSFER VELOCITIES OF CUTTER AND PINION IN CUTTER SYSTEM 
CALL EQVEC(VT1,VT11) 

CALL EQVEC ( VT2 , VT12 ) 

C... RELATIVE VELOCITIS OF CUTTER WRT. PINION IN CUTTER SYSTEM 
CALL ADDVEC(W12,W1,W2,-1.D0) 

CALL ADDVEC(V12 , VT1, VT2 , -1.D0) 

KF=KFP 

KH=KHP 

C... PRINCIPAL CURVATUTES AND DIRECTIONS OF PINION 
CALL CURVT(l) 

C. . . PRINCIPAL DIRECTIONS OF PINION IN PINION SYSTEM 
CALL EQVEC (EFF,ES) 
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CALL EQVEC ( EHH , EQ ) 

C. . . ROTATING VELOCITY OF CUTTER IN CUTTER SYSTEM 
W1(1)=0.D0 
W1(2)=0.D0 
W1(3)=0.D0 

C... ROTATING VELOCITY OF GEAR IN CUTTER SYSTEM 
W2(1)=O.DO 
W2 (2 ) =0 . DO 
W2(3)=-DGPHI2 

C... NORMAL OF CUUTER IN CUTTER SYSTEM 
EX( 1) =XNP2 
EX ( 2 ) =YNP2 
EX(3)=ZNP2 

C... TRANSFER VELOCITIES OF CUTTER AND PINION 
CALL EQVEC (VT1,VT21) 

CALL EQVEC (VT2 , VT22) 

CALL EQVEC (EF,EF2) 

CALL EQVEC (EH ,EH2) 

KF=KF2 

KH=KH2 

C... RELATIVE VELOCITIS OF CUTTER WRT. GEAR IN CUTTER SYSTEM 
CALL ADDVEC (W12 , W1 , W2 , -1 .DO) 

CALL ADDVEC (V12 , VT1, VT2 , -1 . DO) 

KF=KFG 

KH=KHG 

C. .. PRINCIPAL CURVATUTES AND DIRECTIONS OF PINION 
CALL CURVT (2) 

C. . . PRINCIPAL DIRECTIONS OF PINION IN PINION SYSTEM 
CALL TRANSM (L1F,LF1) 

CALL MAVEC (VI , LH2 , ES } 

CALL MAVEC (V2,LFH, VI) 

CALL MAVEC (ES,L1F,V2) 

CALL MAVEC (VI, LH2,EQ) 

CALL MAVEC (V2,LFH, VI) 

CALL MAVEC (EQ,L1F,V2) 

C... CONTACT ELLIPSES 
CALL ELLIP 
ELAL ( I ) =A 

C... AXIS OF PINION ELLIPSE ON TANGENT PLANE 
CALL EQVEC ( VI, SI) 

CALL DOTVEC (V1N , VI , VI) 

ELI (1 , 1) *=DSQRT (VI (1) **2+Vl (2) **2/VlN) 

ELI (2, I) -VI (3) /DSQRT(VIN) 

C... AXIS OF GEAR ELLIPSE ON TANGENT PLANE 
CALL TRANSM (LHF, LFH) 

CALL TRANSM (L2H,LH2) 

CALL MAVEC (V2,LF1, VI) 

CALL MAVEC (VI, LHF, V2) 

CALL MAVEC (V2,L2H, VI) 

CALL EQVEC ( VI, V2) 

CALL DOTVEC (V1N, VI, VI) 

EL2 ( 1 , 1) =-DSQRT ( VI ( 1) **2+Vl (2 ) **2 /V1N) 

EL2 (2,1) =V1 ( 3 ) /DSQRT (V1N) 

AG1=DATAN (ELI ( 2 , 1) /ELI ( 1 , 1 ) ) /All 
AG2»DATAN(EL2 (2,1) /EL2 (1,1)) /All 
DPHI (1, I) =PHI1 
DPHI (2 , I) =PHI2-TN1/TN2*PHI1 
DDR1=DSQRT (RG1 (1) **2+RGl (2 ) **2 ) 

DDR2=DSQRT (RG2 (1) **2+RG2 (2 ) **2 ) -210 . DO 
KM=4 
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0U1=FL0AT(I-1) /FLOAT (KM) 

OUP=AINT(OUl) 

IF(OUl.EQ.OUP) THEN 

WRITE (85, *) DDR1,RG1(3) 

WRITE (86 , *) DDR2 , RG2 ( 3 ) 

C WRITE (87 , *) ELAL(I) , ELI (1, 1) , ELI (2 , 1) , EL2 (1,1) ,EL2 (2,1) 

END IF 

1010 CONTINUE 

1011 CONTINUE 
AD=— 10000 .DO 

DO 1090 I=1,NN,KM 

IF(DPHI(2,I) . NE.O.DO) THEN 
IF(DPHI(2,I) .GT.AD) THEN 
AD=DPHI(2,I) 

ENDIF 
END IF 

1090 CONTINUE 

SS=3 60 . DO /TNI 
DO 1020 J=l, 4 

DO 1030 1=1 , NN , KM 

BD=(DPHI (2,1) —AD) /A11*3600.D0 
DDPHI ( 1 , I , J) =DPHI (1,1) / All+SS* ( J-l ) 

DDPHI ( 2 , I , J) =BD 
1030 CONTINUE 
1020 CONTINUE 
KM=4 

DO 1040 J=l, 4 

DO 1050 1=1 , NN , KM 

IF(DDPHI(1,I,J) .NE.O.DO) THEN 

WRITE(90 , *) DDPHI ( 1, I , J) , DDPHI (2 , I , J) 

ENDIF 

1050 CONTINUE 
1040 CONTINUE 
KM=4 

DO 1060 1=1 , NN , KM 

WRITE (87,*) ELAL(I) , ELI (1,1) ,EL1(2 , 1) , EL2 (1, 1) , EL2 (2 , 1) 
1060 CONTINUE 

WRITE (6,*) '***** PROGRAM FINISHED ****** ' 

STOP 

END 

C 

C ...THE SUBROUTINE IS FOR TCA 
C. 

SUBROUTINE FCNG (X , F , N) 

IMPLICIT REAL*8 ( A-H , O-Z) 

REAL* 8 X(N) , F (N) 

REAL* 8 LF1 (3,3) , LH2 (3,3) , LFH (3,3) 

REAL* 8 RIF (3) ,R2F(3) ,N1F(3) ,N2F(3) ,CO(3) 

REAL* 8 AVC1 ( 3 ) , VTR1 ( 3 ) , AVC2 ( 3 ) , VTR2 ( 3 ) 

REAL* 8 UI ( 3 ) , UJ ( 3 ) , UK ( 3 ) , VI ( 3 ) , V2 ( 3 ) , V3 ( 3 ) ,V4(3) ,V5(3) 

REAL* 8 UP , LAMDP , LAMDB 

REAL* 8 KF2 , KH2 

REAL* 8 KFP , KHP , KFG , KHG 

COMMON /AXIS/ UI , U J , UK 

COMMON /SGI/ RG1,NG1,AP,DP 

COMMON /SG2 / RG2 , NG2 , AA 

COMMON /ATT/ PHI1,PHI2,R1F,R2F,N1F,N2F,LFH,CC,DGAM,DPHI2 
COMMON /AST/ ICONT 
COMMON /MVT/ LF1 , LH2 
REAL*8 KSIG1 , KSIG2 , KFF , KHH 
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REAL*8 EFN ( 3 ) , EHN ( 3 ) ,W1VT2(3) ,WV12(3) ,W2VT1(3) , KF , KH , KS , KQ 

REAL*8 KM2 , KT2 

COMMON /A300/ ES(3),EQ(3) 

COMMON /A310/ KFF,KHH 
COMMON /A340/ EFF(3) ,EHH(3) 

COMMON /A340/ EF1 ( 3) , EH1 (3 ) 

COMMON /A360/ A, B, SI (3) , FI (3) 

COMMON/A200/ Wl(3) ,W2 (3) ,W12(3) ,VT1(3) ,VT2(3) ,V12(3) 

COMMON /A210/ EX(3) ,EF(3) ,EH(3) 

COMMON /A212 / EF2 (3) ,EH2 (3) ,KF2,KH2 
COMMON /A220/ KF, KH 
COMMON /A230/ ET(3) f EM(3) 

COMMON /A380/ KS,KQ 

COMMON /A400/ VT11(3) ,VT12 (3) , VT21(3) , VT22 (3) 

COMMON /A401/ KHP , KFP , KHG , KFG 
REAL*8 RG1(3) ,NG1(3) 

REAL*8 RG2 ( 3 ) , NG2 ( 3 ) 

REAL* 8 LAMDP2 , LAMDB2 
REAL *8 LAMDP1 , LAMDB1 

COMMON /DATT2/ TN2 , PN2 , PSIN2 , BETAP2 , ADG2 , DEG2 , LAMDP2 , 

& UP2 , FW2 , RPT2 , RBT2 , RAT2 , RDT2 , PSIT2 , LAMDB2 
COMMON /B2 /XNP1 , YNP1 , ZNP1 , XNP2 , YNP2 , ZNP2 
COMMON /B4/X1, Y1 , Z1 , XN1 , YN1 , ZN1 
COMMON /B5/X2 , Y2 , Z2 , XN2 , YN2 , ZN2 
COMMON /B6/THET2P, DGPHI2 

COMMON /Wl/ ETAW1 , UPP1 , SPP1 , ETAW2 , UPP2 , SPP2 
COMMON /NET / RR,DD 

COMMON /DATT/ TNI , PN1 , PSIN1 , BETAP1 , ADG1 , DEG1 , LAMDP1 , 

& FW1 , RPT1 , RBT1 , RATI , RDT1 , PSIT1 , LAMDB1 
UPP1=X(1) 

ETAW1=X (2 ) 

UPP2=X ( 3 ) 

ETAW2=X (4 ) 

PHI2«X(5) 

IF(ICONT.EQ.l) THEN 
PHI1=X(6) 

ENDIF 

RPT=RPT1 

CNST=DARCOS(-1.0D00) / 18 0.0 
PI=DARCOS (-1 . ODOO) 

PI2=2.D0*PI 
PSIT= PSIT1 
PSIN= PS INI 
BETAP= BET API 
DSIN1=DSIN (ETAW1) 

DCOSl=DCOS (ETAW1) 

ANF l=DARCOS ( RBT1 / RPT1 ) 

CINV1=DTAN (ANF1) -ANF1 
ANG1=PI / 2 . D0/TN1 
AM=PI/PN1/4.D0 
SA=AM*DCOS (BETAP) 

SS1= RPT1 * ( ETAW1 ) 

DSS1= RPT1 

C. . . EQUATION OF MESHING OF PINION 
FF1=SPP1*DSIN (BETAP) 

FF2=(- (UPP1+DP) —2 .D0*AP**2*UPP1**3 ) 

FF3=DSIN (PSIN) +2 . DO*AP*UPP1*DCOS (PSIN) 

F5=FF2/FF3*DCOS (BETAP) -SA-RPT1*ETAW1 
SPP1= F5/DSIN (BETAP) 

C... SURFACE OF PINION RACK CUTTER 
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XP= (UPP1+DP) *DCOS (PSIN) 

YP=( (UPP1+DP) *DSIN (PSIN) +AM) *DCOS (BETAP) + 

& AP*UPP1**2*DC0S (PSIN) *DCOS (BETAP) +SPP1*DSIN (BETAP) 

ZP=- ( (UPP1+DP) *DSIN (PSIN) +AM) *DSIN (BETAP) - 
& AP*UPP1**2*DC0S (PSIN) *DSIN (BETAP) +SPP1*DC0S (BETAP) 

Xl*= DCOS (ETAW1) *XP+DSIN (ETAW1) *YP+RPTl*DCOS (ETAW1) 

&+SSl*DSIN ( ETAW1 ) 

Y1=-DSIN (ETAW1) *XP+DCOS (ETAW1) *YP-RPT1*DSIN (ETAW1) 

&+SSl*DCOS (ETAW1) 

Z1=ZP 

NORMAL OF RACK CUTTER 
PNN=DSQRT ( 1 . D0+ ( 2 . D0*AP*UPP1 ) ** 2 ) 

XNPP=-2 . D0*AP*UPP1/PNN 
YNPP=1 • DO/PNN 
ZNPP=0 . DO 

XNP1= (DCOS (PSIN) *XNPP-DSIN (PSIN) *YNPP) 

YNP1= (DSIN (PSIN) *DCOS (BETAP) *XNPP+DCOS (PSIN) *DCOS (BETAP) *YNPP) 
ZNP1=(-DSIN(PSIN) *DSIN (BETAP) *XNPP-DCOS (PSIN) *DSIN (BETAP) *YNPP) 
NORMAL OF PINION IN SI 
XN1= DCOS (ETAW1) *XNP1+DSIN (ETAW1) *YNP1 
YN1=-DSIN (ETAW1) *XNPl+DCOS (ETAW1) *YNP1 
ZN1= ZNP1 

DX1=-DSIN (ETAW1) *XP+DCOS (ETAW1) *YP+SSl*DCOS (ETAW1) 

DY1— DCOS (ETAW1) *XP-DSIN (ETAW1) *YP-SS1*DSIN(ETAW1) 

DZ1— O.DO 
RG1 (1)=X1 
RG1(2)«Y1 
RG1 (3) =Z1 
NG1(1)=XN1 
NG1 (2) =YN1 
NG1(3)=ZN1 
KHP=0 . DO 

KFP=2 . DO*AP/ (DSQRT ( ( 1 . D0+ ( 2 . D0*AP*UPP1) **2 ) ) **3 ) 

EFl(l)=DCOS (PSIN) 

EF1 (2 ) =DSIN (PSIN) *DCOS (BETAP) 

EF1 ( 3 ) =-DSIN (PSIN) *DSIN (BETAP) 

EH1(1)=0.D0 

EH1 (2) = DSIN (BETAP) 

EH1 (3 ) — DCOS (BETAP) 

EF( 1) = DCOS (ETAW1) *EH1 (1) +DSIN (ETAW1) *EH1 (2) 

EF ( 2 ) =-DSIN ( ETAW1 ) * EH1 ( 1 ) +DCOS ( ETAW1 ) *EH1 ( 2 ) 

EF ( 3 ) = EH1 (3) 

EH(1)= DCOS (ETAW1) *EF1 (1) +DSIN(ETAW1) *EF1(2) 

EH ( 2 ) =-DSIN ( ETAW1 ) *EF1 ( 1 ) +DCOS ( ETAW1 ) *EF1 ( 2 ) 

EH (3 ) = EF1 (3) 

VELOCITY OF RACK-CUTTER IN SP 
VT11(1)=0.D0 
VT11 (2 ) =RPT1 
VT11 (3 ) =0.D0 

VELOCITY OF PINION IN SP 
VT12 (1)=-(YP+SS1) 

VT12 (2) * XP+RPT1 

VT12 (3) =0.D0 

RPT=RPT2 

PSIT= PSIT2 

PSIN= PSIN2 

BETAP= BETAP 2 

SS2= RPT2 * ( ETAW2 /TN2 *TN 1 ) 

DSS2= RPT2/TN2*TN1 
TTD=AA* ( ETAW2 -THET2P ) **2 
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DTP=2 . DO*AA* ( ETAW2 -THET2P ) 

GPHI2= ETAW2/TN2*TN1-TTD 
DGPHI2= l.D0/TN2*TNl-DTP 
C... EQUATION OF MESHING OF GEAR 

GG1=SPP2*DSIN(PSIN2) *DSIN(BETAP2) 

GG1-DSIN (PSIN2 ) *DSIN (BETAP2 ) 

GG2=(UPP2+DP) *DCOS (BETAP2) 

GG3= RPT2 *TNl/TN2*ETAW2*DSIN (PSIN2) 

GG4=SA*DSIN (PSIN2 ) 

GG5= 2 . D0*AA*TN2* (ETAW2-THET2P) *RPT2*DCOS (PSIN2) *DCOS (BETAP2 ) 
GG6= TN1-2 .D0*AA*TN2* (ETAW2-THET2P) 

SPP2=«* ( GG24GG3 +GG4 +GG5 / GG6 ) /GG1 
C. . . SURFACE OF GEAR RACK CUTTER 
XP= (UPP2+DP) *DCOS (PS IN) 

YP= ( (UPP2+DP) *DSIN (PSIN) +AM) *DCOS (BETAP) +SPP2*DSIN (BETAP) 

ZP=- ( (UPP2+DP) *DSIN (PSIN) +AM) *DSIN (BETAP) +SPP2*DCOS (BETAP) 
X2=-DCOS (GPHI2) *XP+DSIN (GPHI2) *YP+RPT2 *DCOS (GPHI2) 

&+SS2*DSIN (GPHI2) 

Y2— “DSIN (GPHI2) *XP-DCOS (GPHI2) *YP+RPT2*DSIN (GPHI2) 

&-SS2*DCOS (GPHI2) 

Z2=ZP 

C. . . NORMAL OF RACK CUTTER IN SW 
XNP2=DSIN (PSIN) 

YNP2=-DCOS (PSIN) *DCOS (BETAP) 

ZNP2= DCOS (PSIN) *DSIN (BETAP) 

XNP2=-DSIN (PSIN) 

YNP2= DCOS (PSIN) *DCOS (BETAP) 

ZNP2=-DCOS (PSIN) *DSIN (BETAP) 

C . . NORMAL OF GEAR IN S2 

XN2=-DCOS (GPHI2) *XNP2+DSIN (GPHI2) *YNP2 
YN2=-DSIN(GPHI2) *XNP2-DCOS (GPHI2) *YNP2 
ZN2= ZNP2 

DX2= (DSIN (GPHI2 ) *XP+DCOS (GPHI2 ) * YP-RPT2 *DSIN (GPHI2 ) + 
&SS2*DCOS (GPHI2 ) ) *DGPHI2+DSS2 *DSIN (GPHI2 ) 

DY2= ( -DCOS (GPHI2) *XP+DSIN (GPHI2) * YP+RPT2 *DCOS (GPHI2) + 
&SS2*DSIN (GPHI2) ) *DGPHI2-DSS2*DC0S (GPHI2) 

DZ2=0.D0 
RG2 (1) =X2 
RG2 (2) =Y2 
RG2 (3) =Z2 
NG2 (1) =XN2 
NG2 ( 2 ) =YN2 
NG2 (3) =ZN2 
KFG=0 . DO 
KHG=0 . DO 

EF2 (1) =DCOS(PSIN) 

EF2 (2) =DSIN(PSIN) *DCOS (BETAP) 

EF2 (3 ) =-DSIN (PSIN) *DSIN (BETAP) 

EH2 (1) =0.D0 

EH2 ( 2 ) = DSIN (BETAP) 

EH2 (3) = DCOS (BETAP) 

C VELOCITY OF RACK-CUTTER IN SG 

VT21 (1) “0 . DO 
VT21 (2 ) = RPT1 
VT21(3)=0.D0 
C VELOCITY OF GEAR 

VT22 (1) — ( YP+SS2 ) *DGPHI2 
VT22 (2) = (XP+RPT2) *DGPHI2 
VT22 (3)=0.D0 
S1=DSIN(PHI1) 
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C1=DC0S (PHI1) 

S2=DSIN(PHI2) 

C2=DC0S (PHI2 } 

LF1 (1,1) =C1 

LF1 (1, 2) =S1 

LF1 ( 2 , 1) =-Sl 

LF1 (2,2) =C1 

LF1 (1, 3 ) =0 .DO 

LF1 (2 , 3 ) =0 .DO 

LF1 (3 , 1) =0.D0 

LF1 (3 , 2) =0.D0 

LF1 ( 3 , 3) =1.D0 

LH2 (1,1) =C2 

LH2 (1,2) =-S2 

LH2(2,1)=S2 

LH2 (2 , 2) =C2 

LH2 (1,3) —0 . DO 

LH2(2,3)=0.D0 

LH2 (3,1) =0.D0 

LH2 (3,2) =0 . DO 

LH2 (3 , 3) =1.D0 

CALL MAVEC (VI , LH2 , NG2 ) 

CALL MAVEC (N2F,LFH, VI) 

CALL MAVEC (V2,LH2,RG2) 

CALL MAVEC (V3,LFH,V2) 

CALL ADDVEC (R2F , V3 , UI , CC) 

CALL MAVEC (N1F,LF1,NG1) 

CALL MAVEC (RIF, LF1,RG1) 

F(1)=(R1F(1)-R2F(1) ) 

F(2)=(R1F(2) -R2F (2) ) 

F(3) =(R1F(3) -R2F (3) ) 

F(4)=N1F(1) -H2F (1) 

F(5)=N1F(2)-N2F(2) 

IF(ICONT.EQ.l) THEN 
F(6)=Z1-DD 
ENDIF 
RETURN 
END 

C. . . 

C. . . FOR PINION DATA ** 

C. • • 

SUBROUTINE DATAT1 
IMPLICIT REAL*8 (A— H,0-Z) 

REAL *8 UP , LAMDP1 , LAMDB1 

COMMON / DATT/ TN1,PN1,PSIN1,BETAP1,ADG1,DEG1,LAMDP1, 

& FW1 , RPT1 , RBT1 , RATI , RDT1 , PSIT1 , LAMDB1 
RPT1=TN1/ (2 .DO*PN1*DCOS (BETAP1) ) 

PSIT1=DATAN (DTAN (PSIN1) /DCOS (BETAP1) ) 

RBT1=RPT1*DC0S (PSIT1) 

RAT 1 =RPT 1+ADG1 
RDT1=RPT1-DEG1 

LAMDB1=DATAN (DTAN (LAMDP1) /DCOS (PSIT1) ) 

BETAB1=DATAN ( DTAN ( BETAP1) *DC0S(PSIT1) ) 

WRITE (55, 110) 

HO FORMAT ( / 2 X, ' ****************************************** f j 
& 2X , ' * DATA OF PINION *• ,/ 

& 2X , ' it****************************************/ / ) 

WRITE (55,120) TN1,PN1, PSIN1 , BETAP1 , ADG1 , DEG1 , LAMDP1 , FW1 , 

& RPT1 , RBT1 , RATI , RDT1 , PSIT1 , LAMDB1 
120 FORMAT (2X, 'GEAR NUMBER OF TEETH TN1=',F14.7/ 
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&2X, 'NORMAL DIAMETRAL PITCH (1/MM) PN1=',F14.7/ 

&2X, 'NORMAL PRESSURE ANGLE (RAD.) PSIN1*' , F14 . 7/ 

&2X, 'LEADING ANGLE OF HELIX' / 

&4X, 'ON PITCH CYLINDER (RAD.) BETAP1=' , F14 .7/ 

&2X, 'ADDENDUM (MM) ADG1=' , F14 .7/ 

&2X, 'DEDENDUM (MM) DEG1=' , F14 .7/ 

&2X, 'HELIX ANGLE ON PITCH CYLINDER (RAD) LAMDP1*' , F14 . 7/ 
&2X, 'FACE WIDTH (MM) FW1=',F14.7/ 

&2X, 'RADIUS OF PITCH CYLINDER (MM) RPT1=' , F14 . 7/ 

&2X, 'RADIUS OF BASE CYLINDER (MM) RBT1=' , F14 . 7/ 

&2X, 'RADIUS OF ADDENDUM CYLINDER (MM) RATI*' , F14 . 7/ 

&2X, 'RADIUS OF DEDENDUM CYLINDER (MM) RDT1=' , F14 . 7/ 

&2X, 'TRANSVERSE PRESSURE ANGLE (RAD> PSIT1*' , F14 . 7/ 

&2X, 'HELIX ANGLE ON BASE CYLINDER (RAD) LAMDB1*' , F14 . 7/ ) 

C. . . 

RETURN 

END 

C. . . 

C ...THE SUBROUTINE IS FOR DATA OF THEORITI CAL GEAR SURFACE 
C. • 

SUBROUTINE DATAT2 
IMPLICIT REAL*8 (A-H,0-Z) 

REAL* 8 UP2 , LAMDF2 , LAMDB2 

COMMON / DATT2 / TN2 , PN2 , PS IN2 , BETAP2 , ADG2 , DEG2 f LAMDP2 , 

& UP2 , FW2 , RPT2 r RBT2 , RAT2 , RDT2 , PSIT2 , LAMDB2 
RPT2=TN2 / ( 2 . D0*PN2 *DCOS (BETAP2 ) ) 

PSIT2=DATAN (DTAN (PSIN2 ) /DCOS (BETAP2 ) ) 

RBT2=RPT2*DC0S (PSIT2) 

RAT2=RPT2+ADG2 
RDT2=RPT2 -DEG2 

LAMDB2 =DATAN (DTAN (LAMDP2 ) /DCOS (PSIT2 ) ) 

BETAB2=DATAN (DTAN (BETAP2 ) *DCOS (PSIT2 ) ) 

DEL2=0. 007 
C WRITE (6, 110) 

WRITE(55 , 110) 

110 FORMAT(/2X , '*****************************************/ f / 

& 2X, ' * DATA OF GEAR 2 *',/ 

& 2X, ' ****************************************** f j) 

WRITE(55, 120) TN2 , PN2 , PSIN2 , BETAP2 , ADG2 , DEG2 , LAMDP2 , 

& FW2 , RPT2 , RBT2 , RAT2 , RDT2 , PSIT2 , LAMDB2 , DEL2 
120 FORMAT (2X, 'GEAR NUMBER OF TEETH TN2=',F14.7/ 

&2X, 'NORMAL DIAMETRAL PITCH (1/MM) PN2=',F14.7/ 

&2X, 'NORMAL PRESSURE ANGLE (RAD.) PSIN2=' , F14 . 7/ 

&2X, 'LEADING ANGLE OF HELIX' / 

&4X, 'ON PITCH CYLINDER (RAD) BETAP2=' , F14 .7/ 

&2X, 'ADDENDUM (MM) ADG2=' , F14 . 7/ 

&2X, 'DEDENDUM (MM) DEG2=' , F14 . 7/ 

&2X, 'HELIX ANGLE ON PITCH CYLINDER ( RAD ) LAMDP2=' , F14 . 7/ 
&2X, 'FACE WIDTH (MM) FW2=' / F14.7/ 

&2X, 'RADIUS OF PITCH CYLINDER (MM) RPT2=' ,F14 . 7/ 

&2X, 'RADIUS OF BASE CYLINDER (MM) RBT2=' , F14 . 7/ 

&2X, 'RADIUS OF ADDENDUM CYLINDER (MM) RAT2=' , F14 . 7/ 

&2X, 'RADIUS OF DEDENDUM CYLINDER (MM) RDT2=' , F14 .7/ 

&2X, 'TRANSVERSE PRESSURE ANGLE (RAD.) PSIT2=' , F14 . 7/ 

&2X, 'HELIX ANGLE ON BASE CYLINDER ( RAD ) LAMDB2=' , F14 . 7 / 

&2X, 'ELASTIC APPROACH (MM) DEL*' , F14 . 7) 

C. • 

RETURN 

END 

C. . . 
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C... ADDITION OF TWO VECTORS 
C. . . 

SUBROUTINE ADDVEC (VA, VB, VC, DD) 
IMPLICIT REAL *8 (A-H, O-Z) 

REAL* 8 VA ( 3 ) , VB ( 3 ) , VC ( 3 ) 

DO 101 1=1,3 

VA ( I ) =VB ( I ) +DD*VC ( I ) 

101 CONTINUE 
RETURN 
END 

C. . . 

C... DOT PRODUCT OF TWO VECTOR 
C. . . 

SUBROUTINE DOTVEC (AA, VA, VB) 
IMPLICIT REAL*8 (A-H, O-Z) 

REAL*8 VA(3) ,VB(3) 

AA=0 . DO 
DO 102 1=1,3 

AA=AA+VA (I) *VB (I) 

102 CONTINUE 
RETURN 
END 

C. . . 

C. . . CROSS PRODUCT OF TWO VECTOR 
C. . . 

SUBROUTINE CROVEC ( VA , VB , VC ) 
IMPLICIT REAL* 8 (A-H, O-Z) 

REAL*8 VA(3) ,VB(3) ,VC(3) 
VA(1)=VB(2) *VC(3) -VB(3) *VC(2) 

VA ( 2 ) =VB ( 3 ) *VC ( 1 ) -VB ( 1 ) *VC ( 3 ) 

VA ( 3 ) =VB ( 1 ) *VC ( 2 ) — VB ( 2 ) *VC(1) 

RETURN 

END 

C. . . 

C. . . PRODUCT OF MATRIX AND A VECTOR 
C. . . 

SUBROUTINE MAVEC ( VA , MC , VB ) 
IMPLICIT REAL*8 (A-H, O-Z) 

REAL*8 MC (3,3) , VA ( 3 ) , VB ( 3 ) 

DO 103 1=1,3 
VA ( I ) =0 . 0 
DO 104 J-1,3 

VA(I) =MC (I , J) *VB ( J) +VA (I) 

104 CONTINUE 

103 CONTINUE 
RETURN 
END 

C. . . 

C. . . PRODUCT OF A VECTOR AND A SCALAR 
C. • 

SUBROUTINE PDSVEC(VA,VB,T) 
IMPLICIT REAL *8 (A-H, O-Z) 

REAL* 8 VA(3) ,VB(3) 

DO 105 1=1,3 
VA(I) =T*VB (I) 

105 CONTINUE 
RETURN 
END 

C. 

C. . . STANDARDIZATION OF A VECTOR 


69 



C • « • 


SUBROUTINE STDVEC (VA, VB) 

IMPLICIT REAL*8 (A-H,0-Z) 

REAL* 8 VA(3) f VB(3) 

CC=O.DO 
DO 106 1=1,3 
CC=CC+VB ( I ) **2 

106 CONTINUE 
CN=DSQRT(CC) 

DO 107 1=1,3 

VA ( I ) =VB ( I ) / CN 

107 CONTINUE 
RETURN 
END 

C. 

cl 1 1 INPUT A VECTOR TO ANOTHER VECTOR 

C. 

SUBROUTINE EQVEC ( VA , VB ) 

IMPLICIT REAL *8 (A-H,0-Z) 

REAL* 8 VA(3) ,VB(3) 

DO 108 1=1,3 
VA ( I ) =VB ( I ) 

108 CONTINUE 
RETURN 
END 

C. 

dll TRIPLE PRODUCT OF THREE VECTORS 
C. 

SUBROUTINE TRIVEC (AA, VA, VB , VC) 

IMPLICIT REAL* 8 (A-H,0-Z) 

REAL*8 VA(3) ,VB(3) ,VC(3) ,V(3) 

CALL CROVEC(V, VB, VC) 

CALL DOTVEC(AA, VA, V) 

RETURN 

END 

c. 

cl.. TRANSFOR MATRIX 
C. . . 

SUBROUTINE TRANSM(AA, BB) 

IMPLICIT REAL* 8 (A-H,0-Z) 

REAL*8 AA (3,3) , BB (3,3) 

AA (1,1) =BB (1,1) 

AA (1,2) =BB (2,1) 

AA (1,3) =BB (3,1) 

AA (2,1) =BB (1,2) 

AA (2,2) =BB (2,2) 

AA (2,3) =BB (3,2) 

AA (3,1) =BB (1,3) 

AA (3,2) =BB (2,3) 

AA (3,3) =BB (3,3) 

RETURN 

END 

C. . . 

C. . . PRODUCT OF TWO MATRICES 
C. 

SUBROUTINE MAPMA(AA, BB,CC) 

IMPLICIT REAL* 8 (A-H,0-Z) 

REAL* 8 AA (3,3) , BB (3,3) , CC (3,3) 

AA( 1, 1) =BB(1 , 1) *CC (1, 1) +BB (1 , 2) *CC (2 , 1) +BB(1, 3) *CC (3 , 1) 
AA(1, 2) =BB (1,1) *CC(1, 2) +BB (1, 2) *CC (2 , 2) +BB (1, 3) *CC (3 ,2) 
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AA(1 , 3 ) =BB (1,1) *CC (1, 3) +BB (1,2) *CC (2 , 3) +BB (1, 3) *CC(3 , 3 ) 
AA(2,1)=BB(2,1)*CC(1,1)+BB(2,2) *CC (2 , 1) +BB (2 , 3) *CC (3 , 1) 

AA (2,2) =BB (2,1) *CC (1, 2 ) +BB (2 , 2) *CC (2 , 2) +BB (2 , 3) *CC(3 , 2) 
AA(2 f 3)=BB(2 f 1) *CC (1 , 3) +BB (2 , 2) *CC< 2 , 3 ) +BB (2 , 3) *CC (3 , 3) 

AA(3 , 1) =BB (3,1) *CC (1, 1) +BB (3 ,2) *CC(2 ,1) +BB (3 , 3) *CC(3 ,1) 

AA (3,2) =BB (3,1) *CC (1,2) +BB (3,2) *CC (2,2) +BB (3,3) *CC (3,2) 

AA(3 , 3) =BB (3,1) *CC( 1, 3) +BB (3 f 2) *CC (2 , 3 ) +BB (3 , 3) *CC(3 , 3) 

RETURN 

END 

C. . 

C. . . THE SUBROUTINE IS FOR MISALIGMENT 
SUBROUTINE MIAL ( K , DGAMM , LFH ) 

IMPLICIT REAL*8(A-H,0-Z) 

REAL* 8 LFH (3,3) 

S3=DSIN (DGAMM) 

C3=DC0S (DGAMM) 

C ** FOR CROSSING ANGLE MISALIGNMENT 
IF(K.EQ.l) THEN 
LFH (1,1) =~1 • DO 
LFH( 1, 2) =0 .DO 
LFH(1, 3) =0 .DO 
LFH(2,1)=O.DO 
LFH (2,2) =-C3 
LFH (2,3) =-S3 
LFH(3 , 1) =O.DO 
LFH (3,2) =-S3 
LFH (3,3) =C3 
ENDIF 

C ** FOR INTERSECTING ANGLE MISALIGNMENT 
IF(K.EQ.2) THEN 
LFH(1, 1)=-C3 
LFH (1, 2) =0 .DO 
LFH(1 , 3 ) =S3 
LFH (2, 1)=0.D0 
LFH(2,2)*-1.D0 
LFH (2 , 3 )=0.D0 
LFH (3,1) =S3 
LFH (3,2) =0 . DO 
LFH (3,3) =C3 
ENDIF 

C 

RETURN 

END 

C 

C COMPUTE THE PRINCIPAL CURVATURES OF MODIFIED cutter SURFACE 

C 

SUBROUTINE CURVT(KK) 

IMPLICIT REAL* 8 (A-H,0-Z) 

REAL*8 EFN(3) ,EHN(3) ,W1VT2(3) ,WV12(3) ,W2VT1(3) ,KF,KH,KS,KQ 
REAL* 8 KM2 , KT2 , KFF , KHH 

COMMON/ A2 00 / Wl(3) ,W2(3) ,W12(3) ,VT1(3) ,VT2(3) ,V12(3) 

COMMON /A210/ EX (3) , EF (3) , EH(3 ) 

COMMON /A220/ KF , KH 
COMMON /A230/ ET(3),EM(3) 

COMMON /A300/ ES(3),EQ(3) 

COMMON /A310/ KFF, KHH 

COMMON /A380/ KS,KQ 

EFN ( 1) = EX ( 2 ) *EF ( 3 ) -EX (3)*EF(2) 

EFN(2)=-(EX(1) *EF (3) -EX (3) *EF (1) ) 

EFN(3)= EX ( 1) *EF (2) -EX (2) *EF ( 1) 
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EHN(1)= EX (2) *EH(3) -EX (3) *EH(2) 

EHN(2)=-(EX(1) *EH(3) -EX (3) *EH(1) ) 

EHN (3 ) = EX (1) *EH (2) -EX(2 ) *EH( 1) 

W1VT2 (1) = W1 ( 2 ) * VT2 ( 3 ) -W1 ( 3 ) *VT2 ( 2 ) 

W1VT2 ( 2 ) =- ( W1 ( 1 ) *VT2 ( 3 ) -W1 ( 3 ) *VT2 ( 1 ) ) 

W1VT2 ( 3 ) = Wl (1) *VT2 (2) -W1 (2) *VT2 (1) 

W2VT1(1)= W2 ( 2 ) * VT1 ( 3 ) -W2 ( 3 ) *VT1 ( 2 ) 

W2VT1 ( 2 ) =- (W2 ( 1 ) *VT1 ( 3 ) -W2 ( 3 ) *VT1(1) ) 

W2VT1(3)= W2 ( 1 ) *VT1 ( 2 ) -W2 ( 2 ) *VT1 ( 1 ) 

WV12 (1) = W12 (2) *V12 (3) -W12 (3) *V12 (2) 

WV12 (2)=-(W12 (1) *V12 (3) -W12 (3) *V12 (1) ) 

WV12 (3) = W12 (1) *V12 (2) -W12 (2) *V12 (1) 

V12F=0. 0 
V12H=0. 0 
WNEF=0 . 0 
WNEH-0 . 0 
VWN= 0.0 
W1TN=0 . 0 
W2TN=0.0 
DO 1 1=1,3 

V12F= V12 (I) *EF(I) +V12F 
V12H= V12 (I) *EH (I) +V12H 
WNEF= W12 (I) *EFN(I) +WNEF 
WNEH= W12 (I) *EHN (I) +WNEH 
VWN = EX ( I ) *WV12 ( I ) +VWN 
W1TN= EX (I) *W1VT2 (I) +W1TN 

1 W2TN= EX(I) *W2VT1 (I) +W2TN 

C... COMPOTE THE CURVATURE OF THE GENERATED SURFACE 
B13=-KF*V12F-WNEF 
B23=-KH*V12H-WNEH 

B33=-KF*V12F**2-KH*V12H**2+VWN-W1TN+W2TN 

B11=B13**2/B33 

B12=B13*B23/B33 

B22=B23**2/B33 

Tl=2.0DO0*B13*B23 

T2=B23**2-B13**2- (KF-KH) *B33 

SIG1F=0 . 5DOO*DATAN (T1/T2 ) 

C PRINCIPAL CURVATURES OF THE GENERATED SURFACE 

IF (DABS ( SIG1F) . LE . 0 . ID-5 ) THEN 

KQ=0 . 5D0* ( (KF+KH) + (B13**2+B23**2) /B33 
& +(B23**2-B13**2-(KF-KH) *B33) / (B33*DCOS (SIG1F) ) ) 

KS=KQ-(B23**2-B13**2- (KF-KH) *B33) / (B33*DCOS (SIG1F) ) 
ELSE 

KQ=0 . 50D00* (KF+KH) +0 . 5D00* (B13**2+B23**2) /B33 
& +B13*B23 / (B33*DSIN(2 . ODOO* (+SIG1F) ) ) 

KS= KQ-2 . ODOO*B13*B23/ (B33*DSIN (2 . ODOO* (+SIG1F) ) ) 
ENDIF 

SIGSF=-SIG1F 

C PRINCIPAL DIRECTIONS OF THE GENERATED SURFACE 

DO 2 1=1,3 

EQ (I ) = DCOS ( SIGSF) *EH ( I ) -DSIN ( SIGSF) *EF (I) 

2 ES ( I ) = DSIN(SIGSF) *EH(I) +DCOS(SIGSF) *EF(I) 

IF (KK .LT. 2) GO TO 100 

Q3=ES ( 1) *ET (1) +ES (2) *ET(2) +ES (3) *ET (3 ) 
Q4=EQ(1)*ET(1)+EQ(2)*ET(2)+EQ(3)*ET(3) 

IF(Q3.EQ.Q4) THEN 
Q34=DATAN (Q4) 

ELSE 

Q34=DATAN2(Q4,Q3) 

ENDIF 
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KT2=KS*DCOS(Q34)**2+KQ*DSIN(Q34)**2 
KM2=KS*DSIN(Q34) **2+KQ*DCOS (Q34) **2 
GO TO 120 
100 KFF=KS 
KHH=KQ 
120 RETURN 
END 
C 

C . . . . COMPUTE THE DIMENSIONS AND DIRECTIONS OF THE CONTACT ELLIPSE 
C 

SUBROUTINE ELLIP 
IMPLICIT REAL*8(A-H,0-Z) 

REAL*8 KSIG1 , KSIG2 , KFF , KHH , KS , KQ 
COMMON /A300/ ES(3),EQ(3) 

COMMON /A380/ KS,KQ 
COMMON /A310/ KFF, KHH 
COMMON /A340/ EFF(3) ,EHH(3) 

COMMON /A360/ A,B, SI (3) , FI (3 ) 

DEL=0. 007D0 
PI=DACOS (-1 . DO) 

S1=ES ( 1) *EFF ( 1) +ES ( 2 ) *EFF ( 2 ) +ES ( 3 ) *EFF (3 ) 

S2=EQ (1) *EFF (1) +EQ(2) *EFF (2 ) +EQ (3) *EFF(3) 

SIGSF=DATAN2 (S2 , S1J 

C COMPUTE TfiE DIMENSIONS OF THE CONTACT ELLIPSE (A & B) 

KSIG1=KFF+KHH 

KSIG2=KS+KQ 

G1=KFF-KHH 

G2=KS-KQ 

A= (KSIG1-KSIG2- <Gl**2-2 . DO*G1*G2*DCOS (2 . D0*SIGSF) +G2**2 ) **0.5) 
& / 4 .DO 

B=(KSIG1-KSIG2+ (Gl**2-2 ,DO*G1*G2*DCOS (2 . D0*SIGSF) +G2**2) **0.5) 
& /4 . DO 

A= (DEL/ABS (A) ) **0 . 5 
B= (DEL/ABS (B) ) **0 . 5 

C COMPUTE THE ANGLE (ALF1) BETWEEN AXES OF ELLIPSE & PRINCIPLE 

S1=G2 *DSIN ( 2 . D0*SIGSF) 

S2=G1-G2*DC0S (2 .D0*SIGSF) 

ALF1=0 . 5D0*DATAN2 (S1,S2) 

S3=DSQRT(Gl**2-2 .DO*G1*G2*DCOS (2 .D0*SIGSF) +G2**2) 

SS2=S1/S3 

SC2=S2/S3 

ALF1=DATAN (SS2 / ( 1 . D0+SC2 ) ) 

ALF1=DABS (ALF1) 

C. . AXES OF THE CONTACT ELLIPSE 
DO 100 1=1,3 

SI { I ) =DSIN ( ALF 1) *EFF(I) +DCOS ( ALF 1 ) * EHH ( I ) 

FI ( I ) =DCOS ( ALF1 ) *EFF (I) -DSIN ( ALF 1 ) * EHH ( I ) 

100 CONTINUE 
RETURN 
END 
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